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We present the study of the semileptonic decay A c — > A£ + u by using the light-cone sum rule 
approach. Distribution amplitudes(DAs) for the A baryon are discussed to the leading order con- 
formal spin, and QCD sum rule estimate for the corresponding parameters is presented. The form 
factors describing the decay are calculated and used to predict the decay width and the decay 
asymmetry parameter a. With the inclusion of twist-3 contributions the calculated decay width 
T = (7.2 ± 2.0) x 10 _14 GeV as well as asymmetry a = —(0.88 ± 0.03) is found in good agreement 
\Q . with the experimental data, while there are appreciable deviations from experiment values when the 

' higher twist contributions are included. 
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<^ I. INTRODUCTION 

a\ : 

The study of flavor changing decays of the c-quark is always an active field in the heavy flavor physics, and the 
most obvious reason lies on that those processes can provide useful information on the various charm related Cabibbo- 
Kobayashi-Maskawa(CKM) matrix elements which are the main ingredients of the standard model(SM). Furthermore, 
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the thorough understanding of the SM in itself needs a comprehension of the flavor changing dynamics. Unfortunately, 
such a comprehension is difficult contemporarily, the fact is that form factors characterizing those processes are not 
perturbative quantities and whose determination must invoke some non-perturbative method. This paper aims to 
give a preliminary determination of the form factors of semileptonic A c — > A£ + i> decay. In the calculation we will 
use the method of QCD sum rules on the light-cone which in the past has been successfully applied to various 
problems in heavy meson physics, see Q for a review. 
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The method of light-cone sum rules (LCSR) is a new development of the standard technique of QCD sum rules a 
la SVZ sum rules [3j, which comes as the remedy for the conventional approach in which vacuum condensates carry 
D . no momentum 4]. The main difference between SVZ sum rule and LCSR is that the short-distance Wilson OPE in 
increasing dimension is replaced by the light-cone expansion in terms of distribution amplitudes of increasing twist, 
m ^ ' originally used in the description of the hard exclusive process Q . In recent years there have been manyapplications 
of LCSR to bary ons. The nucleon electromagnetic form factors were studied for the first time in [y, [?J and, more 
recently, in 0, Q for a further consideration. Several nucleon related processes gave fruitful results within LCSR, the 
weak decay Af, — » p£v£ was considered in llfl in both full QCD and HQET LCSR. The generalization to the A^A 
transition form factor was worked out in . 

In this paper we will adopt the LCSR approach to study the exclusive semileptonic decay A c — * A£ + i>. This 
transition had been studied in the literature by several authors, employing flavor symmetry or quark model or both 
in Refs. [13. \l^L H3. Eil lu ^ . There are also QCD sum rule description of the form factors 01, upon which the total 
decay rate are obtained. 

The paper is organized as follows: The relevant A baryon DAs are first discussed in Sec. |H] Following that Sec. 
IIIII is devoted to the LCSRs for the semileptonic A c — > M + v decay form factors. The numerical analysis and our 
conclusion are presented in Sec. IIVI 



II. THE A BARYON DISTRIBUTION AMPLITUDES 



Our discussion in this section for the A baryon DAs parallels with that for the nucleon [18|, so we only list the 
results following from that procedure and for details it is recommended to consult the original paper. The A baryon 



2 



DAs are defined through the matrix element 

A{Q\e ijk v}J ai x)dUa 2 x)s k {a 3 x)\P) = {Ax + 



2 M 2 



AfWisCUpk-, + A 2 M{p l5 C) af} {tA) n 



+ A 3 M{^ 5 C) af3 {^A) 7 + A 4 M 2 ^ l5 C) al3 A 7 + A 5 M 2 { W C) af3 {ia»- x v A) n + A 6 M 3 {^C) aP {^A) 7 , (1) 

where the A generically designates the spinor for the A baryon with momentum P. Only axial-vector DAs are 
presented here, for those with other Lorentz structures do not contribute in the final sum rules. The twist of those 
calligraphic DAs is indefinite, but they can be related to the ones with definite twist as 



-Ax + A 2 - A 3 , 
-2Ax - A 3 - A 4 



■2A, 



Ax=Ax, 2P-xA 2 
2A 3 = A 3 , 4P-xA 4 
AP ■ xA 5 = A 3 - A A , {2P-x) 2 A e = A 1 -A 2 + A 3 + A i -A b + A 6 . (2) 

The twist of A, is given in Tab. Each distribution amplitudes F = Ai can be represented as Fourier integral over 



twist-3 


twist-4 


twist- 5 


twist- 6 


Aj 


A 2 ,A 3 


A 4 ,A 5 


A 6 



TABLE I: The twist for A t . 

the longitudinal momentum fractions x\, x 2 , x 3 carried by the quarks inside the baryon with T^Xi = 1, 



F{a t P -x) = J Vxe- ip - xTliXiai F{ Xi ) . 



The integration measure is defined as 



Vx 



dxxdx 2 dx 3 6(xx + x 2 + x 3 — 1). 



As elucidated in [l9| , those distribution amplitudes are scale dependent and can be expanded into orthogonal funct 
with increasing conformal spin. To the leading conformal spin, or s-wave, accuracy the expansion reads |l8| 

Ax{xi,fj,) = -120xxx 2 x 3 (/)l{n), 

A 2 (x t ,fj,) = - 24a; ix 2 04 (m), 

M{Xi,n) = - I2x 3 {\ ~ x 3 )ijjl{ij), 

A 4 (x t ,n) = -3{1-x 3 )i/)°{ij,), 

A 5 (xi,{i) = - 6x 3 05 (/") 

A 6 { Xi ,n) = -20° On), 



(3) 



(4) 



ions 



(5) 



where the constraint A{xx, x 2l x 3 ) — A{x 2 , x\, x 3 ) arising from the condition that the A baryon has isospin has been 
used in the derivation. All the 6 parameters involved in Eq. J2Jl can be expressed in terms of 2 independent matrix 
elements of local operators. Those parameters are 



-2( A i-/a) 



The normalization of A\ at the origin defines the nucleon coupling constant /a, 

(0 | e ijfe [ U l (0)C 75 ^(0)]^ fc (0) I P) = Iaz ■ PtA{P). 

The remaining parameter Ai is defined by the matrix element 

(0 | t ijk [u\Q)C 1 ^dP{Q)]^s k {Q) | P) = \xMA(P). 



(6) 
(7) 



The twist of the order 0{x 2 ) correction starts from twist five, which is apparent in the definition Due to the 
numerically small contribution it gives in the previous applications of light-cone QCD sum rules [ldj . we do not 
consider it in the following analysis. 
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III. A c -> M + v DECAY FORM FACTORS FROM LIGHT-CONE SUM RULES 



Complying with the standard philosophy in the sum rule analysis, we consider the following correlation function 

z v T v {P,q) = iz v [ d 4 xe^ x (0 \ T{j Ac (())>(*)} I P), (8) 



where j'a c = £ijk{u l Cj5j!;d : >)j!;c k is the current interpolating the A c baryon state, j v = c"7„(l — 7s)s is the weak current, 
C is the charge conjugation matrix, and i, j, fc denote the color indices. The auxiliary light-cone vector z is introduced 
to project out the main contribution on the light-cone. The interpolating current used here is not the unique one, as 
exemplified in case studies 20] for the applications of QCD sum rules, and there can be other choices. The coupling 
constant of the interpolating current to the vacuum can thus be defined as 

(0 | jh c | A C (P')} = f\ c z ■ P'^A C (P'), (9) 

where A C (P') and P' is the A c baryon spinor and four- momentum, respectively. Form factors are given in the usual 
way 



(A c (P-q) \j„\ A(P)) = A c (P-q) 



hlv - i 



h 



<72 

Si7 " +, Ma" <7 "'' ' 



a(p), 



(10) 



in which Ma c is the A c mass, A(P) denotes the A spinor, satisfying PA(P) = MA(P), where M is the A mass and 
P its four-momentum. Those form factors give no contribution in the case of massless final leptons are omitted here. 
Giving those definitions, the hadronic representation of the correlation function JSJl can be written as 

2/a c 



z"T„ 



Ml 



p/2 



[z-P' 



flt + fi- 



.91/ 



52" 



75 



A(P) 



(11) 



'M Ac V M Ac 

where P' = P — q and the dots stand for the higher resonances and continuum. While on the theoretical side, at 
large Euclidean momenta P' 2 and q 2 the correlation function © can be calculated perturbatively and the result, in 
the leading order of a s , is 

d A k z ■ k 



z v T„ 



-2(CVW(i - 7 5 ) M / <Tx 



(2tt) 4 k 2 - m 2 

where m c is the c-quark mass. Substituting (yi into Eq. I|12l) we obtain 



l{k+qyx (0 | efcX(<H(0>*OiO I P) , (12) 



z u T v 



-2(z-P) 2 



+2(z-P) 2 



dx , x 4^H + M 2 



k 2 



dxa 



M / dxq 



X3B3 (x 3 ) 



(fc 2 

2M 3 / dx* 



\Bi(x 3 ) 

-,2\2 



2M 4 / dx 



xlB 2 (x 3 ) 



xlB 2 (x 3 ) 
(fc 2 — m 2 ) 3 

^(1 + 75 )A(P) + - 



^(1- 75 )A(P) 



(13) 



(fc 2 — m 2 ) 2 J (fc 2 — m 2 ) 3 _ 

where fc = X3P — g and the ellipses stand for contributions that arc nonleading in the infinite momentum frame 
kinematics P — > 00, g ~ const., z ~ 1/P- The functions Pi are defined by 

B = j dx x Ai{xi, 1 - x\ - x 3 ,x 3 ), 
Jo 

Bi = -2Ai + A 2 - A 3 - A 4 + A 5) 
B 2 = At - A 2 + A 3 + 1 4 - A 5 + A 6 , 
B 3 = -Ax+A 2 -A 3 . 

The DAs with tildes are defined via integration as follows 



(14) 



A(x 3 ) = 
i(x 3 ) = 



dx'v 



dx'o 



l—Xn 



dxiA(xi, 1 - xi - x 3 ,x' 3 ), 



dx\A( X \ 1 1 — x\ — x 3 , x 3 ). 



(15) 
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The origin of those functions can be traced back to the partial integration adopted to eliminate the factor X/P ■ x 
which appears in the distribution amplitudes. When the next-to-leading order conformal expansion is considered, the 
surface terms completely sum to zero. The term B corresponds to the leading twist contribution. The form factors 
f 2 and g 2 in (|13|l are characterized by the higher twist contributions. 

Equating (|llfl and (|13fl . adopting the quark-hadron duality assumption and employing a Borel improvement of P' 2 
on both sides lead us to the desired sum rules for the form factors /i and fe, 



and 



where 



M 2 x 2 e- S °/ M B 



dx2 e 



-s'/Ml 



Bo 



M 2 



Ml 



M 2 



Try [-Bi{x3) + —j- B 2 (x 3 ) 



M 2 B 



Q 2 



x 2 M 2 



M 2 



Bi(x ) - j^ t xqB 2 {x ) 



M 2 e -s /Ml x 2 



m 2 + Q 2 + x 2 M 2 dx \m 2 + Q 2 + x 2 M 2 



M 2 x 2 B 2 {x ) 



faj2 c -Ml /Ml 



1 



M Aa M 
B 3 (x ) 



1 dx 



3 e -s'/Ml 



M 2 , . 
j^x B 2 (x ) ) + 



B Jx x 3 

M 2 e -s a /Ml x 2 



M 2 

B 3 (x 3 ) - jf 2 -B 2 {x 3 ) 



x e 



-so/Ml 



Q 2 



x 2 M 2 



x B 2 (x ) 



+ x\M 2 dx \m 2 + Q 2 + x\M 2 



s' = (1 - x)M 2 



m 2 c + (l-x)Q 2 



and xq is the positive solution of the quadratic equation for s' 



x 

so- 



(16) 



(17) 



(18) 



2M 2 x = y/(Q 2 + so - M 2 ) 2 + 4M 2 (Q 2 + m 2 ) - (Q 2 + s - M 2 ). 



(19) 



As the sum rules for the form factors g\ and g 2 are identical with those for the /i and f 2 , f\ — g\ and f 2 ~ g 2 , we 
will only discuss the results for f\ and f 2 in the following sections. 



IV. NUMERICAL ANALYSIS AND THE CONCLUSION 



A. Values for /a and Ai 



Parameters /a, Ai appear in the conformal expansion of the DAs, so we have to determine their values before 
proceeding to analyze the LCSRs. According to their definitions, we consider correlation functions 

n(q 2 ) = % J d*xe iqx (0 | T{J t (x)Jj{Q)} | 0). (20) 

where J 4 's are given in © and (JJJ. Following the standard QCD sum rule philosophy, an estimate for /a,Ai and their 
relative sign is straightforward 

2 rso 



4 f 2-M 2 /M 2 B _ 



s(l - xfe- s / M *d S - - f S " x(l - x)(l - 2x)e- s ' M *—, 
3 ./„,2 s 



(21) 



4(2ir) 4 \ 2 M 2 e- M2 / M i 



[(1 - x 2 )(l -8x + x 2 ) - 12x 2 In x] e~ s/M ^ds 



4: PCI - x) 2 e- s ^ds - \ a 2 e- M ^ M -, 

12 J ml 3 



(22) 



4(2^) 4 / A A 1 Afe- M2 / M « = ^ f s K 1 " a; )( 3 + lSx ~ 5 ^ + + 12x111 x ) e~'^da 

6 Jml 



12 



(1-x) 



1, w 2 X 
1+3(1- x)(5--) 



-s/M\ 



ds 



(23) 
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where x — m 2 s /s and m s is the s-quark mass. The sum rule for the f\ has been obtained before [10J where the 
corresponding heavy quark limit is also derived. At the working window s Q ~ 1.6 2 GeV 2 and 1 < Mj, < 2 GeV 2 the 
numerical value for the coupling constant reads 

/a = 6.1xl(T 3 GeV 2 , Ai = -1.2xlO~ 2 GeV 2 . (24) 

The relative sign of Xi/f\ is obtained from the sum rule (|23[l . In the numeric analysis, the standard values a = 
-(2n) 2 (qq) = 0.55GeV 3 , b = (2ir) 2 {a s G 2 /it) = 0.47 GeV 4 and m s = 0.15GeV are adopted. It should be noted that 
our value for /a here does coincide with that obtained in [2l| . 

The sum rule for the coupling constant of A c to vacuum is similar to that for /a, where the simple substitution 
m s — » m c should be made, and the numerical value is f\ o = (6.4 ± 0.7) x 10~ 3 GeV 2 , taken from the interval 1 < 
M| < 2GeV 2 with s ~ lOGeV 2 . 



B. Analysis of the LCSRs 

In the numerical analysis for the form factors, the charm quark mass is taken to be ro r = 1.41 GeV 22], and the 
other relevant parameters, A c and A baryon masses and the value of \V C s\, can be found in 23]. We start with analysis 
of the twist-3 sum rules, in which only twist-3 DA is kept. Substitute the above given parameters into the LCSRs 
and vary the continuum threshold within the range so = 7 — 9 GeV 2 , we find there exist an acceptable stability in 
the range Mg = 5 — 7 GeV 2 for the Borel parameter. The Mg and the q 2 dependence for the corresponding form 
factors are shown in Figs. 1 and 2, respectively. Also given in Fig. 2 are two leading twist results, corresponding to 
only retain Bq in the sum rules. It is apparent that in that approximation, only fa and fa survive. Apart from the 
leading twist DA we discuss in (JSJ, there still exists another form from Chernyak, Ogloblin and Zhitnitskv[2l| : 

A i ° Z {xi) = -21(^ as [0.52(a; 2 + x 2 2 ) + 0.34x 2 - 2.05x1X2 - 0.48x 3 (xi + x 2 )], (25) 

where (p as = 120xia;2a;3 is the asymptotic DA. The corresponding result is also illustrated for contrast. 

For the up to twist-6 sum rules, the stability is agreeable within the range sq = 9 — 11 GeV 2 and M§ = 7 — 9 GeV 2 . 
In that working region, the twist-3 contribution to f\ is the dominant one, amounting to over 90%. However, the case 
is different for fa: the main contribution comes from the twist-4 DAs and its magnitude is approximately ~ 1.5 of the 
twist-3 one in the whole dynamical region, but with a different sign. On account of the relatively small momentum 
transfer, the asymptotic behavior of DAs may not be fulfilled and we need incorporate higher conformal spin in 
the expansion for them. Furthermore, QCD sum rule tends to overestimate the higher conformal spin expansion 




5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 

Mb (GeV 2 ) Mb (GeV 2 ) 

FIG. 1: The dependence on M% of the LCSRs for the form factors /i and /2 at q 2 = 0. The continuum threshold is so = 8GeV 2 
for the twist-3 result and so = lOGeV 2 for the up to twist-6 one. 
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FIG. 2: The dependence on q of the LCSRs for the form factors /i and /2. The "COZ" denotes the result obtained from the 
COZ DAs. The continuum threshold and the Borel parameter are s = 8GeV 2 , M% = 6GeV 2 and so = lOGeV 2 , M% = 8GeV 2 
for the twist-3 and the up to twist-6 results. For the two leading twist results, they are chosen to be the same with those for 
the twist-3 sum rule. 



parameters j^ij, and the corresponding parameter will enter in the coefficients of the higher conformal spin expansion, 
which is well known as the Wandzura-Wilczek type contribution. So at the current stage, the strategy to stay with 
twist-3 sum rules seems to be a good choice. The M B and the q 2 dependence for the corresponding form factors are 
also shown in Figs. 1 and 2. 

Both the form factors in the whole kinematical region, < q 2 < {M\ c — M) 2 , can be fitted well by the dipole 
formula 



fM 



a 2 (q 2 /M 2 ) 2 + ai q 2 /M 2 +l' 



(26) 



and the uncertainties are negligible. Below in Table |n] we give those coefficients for two sets of parameters: M B 
6 GeV 2 , s = 8 GeV 2 for the twist-3 sum rule and Ml = 8 GeV 2 , s = 10 GeV 2 for the up to twist-6 one. 





twist-3 


up to twist-6 


h 
h 


a 2 a\ fi(0) 
1.595 -2.203 0.449 
2.992 -3.329 0.193 


a-z ai fi(0) 
0.993 -1.712 0.392 
0.238 -1.339 -0.083 



TABLE II: The dipole fit for the form factors fi and f 2 with M% = 6 GeV 2 , 
Ml =8 GeV 2 , so = 10 GeV 2 for the up to twist-6 sum rules, (HD an d (Q. 



so 



8 GeV for the twist-3 results and 



Using the obtained form factors, we can calculate the differential decay rate and the total decay width for the decay 
A c — > hl + v. The differential decay rate is shown in Fig. 3. If only the twist-3 amplitude A\ is retained, we have for 
the total decay width F = (7.2 ± 2.0) x 10~ 14 GeV, which agrees well with the data given by the Particle Data Group 
[3]. This result is also in agreement with the QCD sum rule predictions made in 0|- As for the decay asymmetry 
parameter a defined in [25j . we obtain a = —0.88 ± 0.03, which corresponds to the ratio at zero momentum transfer 
/2(0)//i(0) = 0.44 ± 0.05. That value lies very close to the recent experimental data from CLEO 26]. Note that the 
errors quoted above reflect the uncertainty due to the Borel parameter Mb and the continuum sq. The uncertainty 
due to the variation of the other QCD parameters is not included, which may reach 5% or more. 

For a comparison, the total decay width computed from the up to twist-6 form factors is T = (6.8 ± 2.0) x 
10 _14 GeV, whose agreement with the experimental value is good, too. However, for the asymmetry parameter a, 
we get a = —(0.54 ± 0.02), which lies above the particle data group's average [23j and still greater than the latest 
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0.2 0.4 0.6 0.8 1 1.2 



q 2 (GeV 2 ) 

FIG. 3: Differential decay rate for A c — > A£ + u. The parameters follow those in Fig. 2. 

experimental measurement 26] . The ratio of the two form factors at zero momentum transfer of that asymmetry is 
/2(0)//i(0) = —0.22 ± 0.03. This phenomenon, i.e., the twist-3 result agrees better with the experiments, may be 
attributed to our incomplete inclusion of the higher conformal spin components in the expansion for the DAs. 

To summarize, we have given a preliminary investigation on the semileptonic decay A c — > A£ + v using the LCSR 
method. Sum rules for the form factors are derived and used to calculate the decay width and the asymmetry 
parameter. The decay width agrees well with the experimental data both for the twist-3 and the up to twist-6 sum 
rules, while the asymmetry's agreement is not so good: the twist-3 DA alone can account for the experimental value 
but the result including higher twist contributions is not so good. This is partly due to the interplay of our incomplete 
inclusion of the higher conformal spin contributions and the QCD sum rule over-estimated value for Ai. 
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